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1.0 TECHNICAL PROGRESS REPORT 


Our research during this quarter was focused on 1) the preparation and 
presentation of technical material for several meetings and docunents; 2) 
the continuation of the geometric analysis of TM film products; and 3) the 
Initial Investigation of the quality of TM data acquired over our forestry 
study area and processed through the new Thematic Mapper Image Processing 
System (TIPS). 

1.1 Preparation and Presentation of Technical Material 

1.1.1 PECORA VIII Symposium 

Based on the invitation to present our TM image quality 
research results at the PECORA VIII Symposium, a paper was pre- 
pared and presented during the symposium (October 5, 1983). 

The paper was compiled from previous reports and will be pub- 
lished in the proceedings of the Symposium during 1984. 

1.1.2 Development of General Electric Brochure Material 

At the request of Mr. David J. Wright, General Electric 
(Space Systems Division) we prepared narrative and illustrative 
material for the new Landsat-4 brochure being prepared by General 
Electric. Our early work illustrating the usefulness of TM Band 5 
for discriminating crop types on single-date Imagery is being high- 
lighted in the agricultural section of the brochure. In addition 
to photographic illustrations, generalized graphs were prepared 
showing the relative spectral response for several cover types for 
each reflective TM band. A short narrative was drafted Indicating 
the improved quality of the TM data and its usefulness in agri- 
cultural surveys. 

1.1.3 Preparation of IEEE Geoscience and Remote Sensing Paper 

At the invitation of Dr. Vincent V. Salomonson, serving as 
Volume Editor, we prepared a manuscript highlighting our results 
to date on the analysis of TM and MSS Image quality focusing on 
one aspect, namely, the spectral variability of selected agri- 
cultural and forest cover types in California. The paper surmna- 
rized our initial work within a small agricultural area in the 
San Joaquin Valley on the Landsat-4 scene acquired December 8, 

1982, and within the NASA-Bucks Lake Forestry Study Site in the 
northern Sierra Nevada on the Landsat-4 scene acquired August 12, 

1983. The early results from the forestry study site are Included 
in this quarterly report In Section 1.3. The special issue of the 
IEEE Transactions on Geoscience and Remote Sensing highlighting 
early evaluations of Landsat-4 TM and MSS data Is scheduled for 
April 1984. 
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1.1.4 Third Landsat-4 Investigator's Workshop 

Material was prepared and presented at the Third Landsat-4 
Investigator's Workshop on December 6-7 at NASA-Goddard Space 
Flight Center. The presentation was designed to highlight 
only the progress made to date since the last Investigator's 
Workshop. Since we had previously presented material at the 
Early Results Symposium, the ASP Annual Convention and the 
PECORA VIII Symposium, only those results recently obtained 
from the forestry study site were presented. Tabulated sum- 
mary statistics and photographic Images of both the TM data 
and color aerial oblique photographs taken simultaneously to 
the Landsat-4 overpass were presented. These data and se- 
lected color Images presented at this workshop are reproduced 
In Section 1.3. 

1.1.5 Presentation to Local Groups and Visitors 

During this reporting period several presentations on our 
recent results were given at meetings of local professional 
organizations and to distinguished visitors of both the De- 
partment of Forestry and Resource Management and the Space 
Sciences Laboratory. Highlights of these presentations were 
the viewing of high quality photographic images created from 
both the TM and MSS digital data by the EROS Data Center, the 
IBM Palo Alto Scientific Center, and NASA-Goddard Space Flight 
Center. A detailed listing of these presentations is provided 
in Section 2.2. 

1.2 Geometric Analysis of TM and MSS Film Products 
1.2.1 Background 

The geometric quality of the TM and MSS film products are 
being evaluated by making selective photo measurements such as 
scale, linear and area determinations; and by measuring the 
coordinates of known features on both the film products and map 
products, and then relating these paired observations using a 
standard linear least squares regression approach. The major 
emphasis of our work is to analyze the film products from 
Landsat-4 that are generally accessible to the user comnunity. 
There were three types of Landsat-4 film products being gen- 
erated at the EROS Data Center (EDC) and the Goddard Space 
Flight Center (GSFC); (1) standard multispectral scanner (MSS) 
film products, (2) "interim" Thematic Mapper (TM) analytical 
film products, and (3) LAS-Scrounge TM "engineering" film 
products. The standard MSS film products are generated tt the 
EDC using the CCT-PM digital data with a Laser Beam Recorder 
(LBR) to produce master film copies for black-and-white and 
color composite reproduction. The "interim" TM analytical 
film products were generated at EDC using the return beam vldl- 
con (RBV) image production system to make the first generation 
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working masters during the Scrounge environment and prior to the 
operational film generation under the TIPS environment. The LAS- 
Scrounge TM "engineering" film products were generated at 6SFC for 
engineering purposes, archiving, and routing with CCT orders for 
LIDQA Investigators. 

The Interim TM film products are being used for our analysis 
because they most closely represented those products that would 
be used by the user community and were the only film products 
available at the time of our Investigation. Originally, these 
products were to be generated by the LAS-Scrounge during the 
pre-TIPS environment, but as the demand for these products In- 
creased, arrangements were made to have the EDC produce the 
products using' the RBV Image production system. Using the CCT- 
PT data, EDC produced the film masters using the LBR and the 
supporting computer system formerly dedicated to RBV film pro- 
duction. In order to adapt the TM data to this system, the re- 
sulting "Interim" TM analytical fljm product represented only a 
sub-area of a full TM scene produced under operational conditions 
(Figure 1). The Interim TM film product represents approximately 
72 percent of the area of a full TM scene. Th<> format Is 5322 x 
5322 pixels, centered on the full scene, with a resulting Image 
of 20.2 x 20.2 cm coverlnq a land surface area of approximately 
2,300,590 ha. (5,684,760 ac.). 

1.2.2 Approach 

A seven step procedure Is being used to evaluate t ■» geometric 
properties of the TM film products. 

(1) Select the image product to be analyzed. A two times 
enlargement of interim Band 7 image T0318-007 (Path 044, Row 033) 
covering the southern Sacramento Valley on 1 February 1983 will 
be selected for this study. This particular scene Is selected 
because (1) It was currently available, (2) with the exception of 
small portions of the image, It was cloud free, (3) this Band 
provided a sharp Image of moderate contrast, and (4) the area 
covered represented a wide range of land use and elevatlonal zones. 

The major limitation of using this image was that due to the low 
sun elevation at the time of Image acquisition (26 degrees) many 
of the steep canyons In the wildland areas In the western and 
northern portions of the scene were deeply shadowed. We felt 
that this limitation was outweighed by the facts that (1) the 
Image In Its two time enlarged format was available, and (2) all 
the small water bodies were at the maximum water levels which 
would allow for better precision In selecting control points. 

(2) Grid the Image Into nine equal area blocks. The TM image 
will be grldded to ensure that the control points will be evenly 
distributed throughout the Image, and that they will represent all 
elevation zones and land use patterns present In the Image (Figure 2). 
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Scene T031 8-007, Path 044, Row 033, 01 February 83 
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Figure 2. Schematic of the grldded TM scene. The scene 
was grldded Into nine equal area blocks to 
ensure that the control points would be evenly 
distributed throughout the Image. For each 
block, 40-50 control points were selected. 
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(3) Select control points. The control points selected 
to represent those natural and cultural features that can be 
located reliably on both the TM Image and on United States 
Geological Survey 7 1/2* quadrangle maps. In the agricultural 
areas, these features will be predominantly field Intersections 
or Irrigation ditches; In the urban areas they will be predom- 
inately major road Intersection or airfield runways; In the 
wildland areas they will be water bodies, stream courses, and 
converging points of ridges and canyons Each control point 
will be pin pricked on thi Image with a corresponding anno- 
tation made on the map sheet with a .30irm pen. Approximately 
40 to 50 control points will be selected for each block on the 
TM Image, and they will be located on 144 map sheets that rep- 
resent the area of the Image. 

(4) Measure Image and map coordinates for control. The 

x and y coordinates of the Image control points will be measured 
to the nearest 0.001 Inch (25.4 x 10"&m.). The corresponding 
UTM east and north map coordinates will be scaled off the map 
sheets to the nearest 10 meters ground distance. 

(5) Check for Image and map coordinate errors. In order 
to check for digitizing and map scaling errors, a first order 
regression between map and Image coordinates will be performed 
using the program developed by Daniel (1971). For those co- 
ordinate pairs for which the residuals are excessive, the digi- 
tized and map coordinates will be verified, and changed where 
appropriate. 

(6) Develop regression between Image and map coordinates. 
Using the odd numbered control points, first and second order 
regressions will be developed between the Image and map co- 
ordinates. 

(7) Evaluate the geometric properties of the image. Using 
the regression developed In (6) the map coordinates of the even 
numbered control points (test points) will be predicted using the 
corresponding digitized Image coordinates. The residuals, sum- 
marized by blocks, will be plotted; and an analysis of variance 
(AN0VA) will be performed on the residuals using the nine blocks 
as treatments. In addition, first and second order regressions 
between the Image and map test points will be developed and the 
resulting coefficients will be compared with those calculated 

In (6). 

1.2.3 Results to Date 

At this time, the first five steps of this procedure have been 
completed. A summary of this work Is given In Table 1. A total 
of 476 control points have been selected. Of these, 238 have been 
designated as "test points". These points are currently being 
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Table 1. 

Summary of 

Control Point Selection for Scene T0318-007 (Path 044, Row 033), 01 February 1983 



BLOCK 

CLOUD 
COVER (£) 

ELEVATION RANGE 
OF POINTS (FEET) 

NUMBER OF 
REGRESSION 
POINTS 

NUMBER OF 
TEST 
POINTS 

GEOGRAPHIC DESCRIPTION 

PREDOMINANT 
LAND USE* 

AA 

10 

14-2080 

27 

27 

Napa Valley; Coast Range 

3, 

1. 4 

AB 

50 

5-384 

17 

17 

Coast Range, Sacramento Valley 

1, 

3, 4 

AC 

5 

0-391 

22 

22 

Sacramento Valley; Sierra Nevada foothills 

1, 

3, 4 

BA 

15 

43-1960 

30 

30 

Coast Range; Sacramento Valley 

3, 

1 

BB 

0 

20-199 

27 

27 

Sacramento Valley 

1 


BC 

0 

20-2000 

27 

27 

Sacramento Valley; Sierra Nevada foothills 

3, 

1, 4 

CA 

0 

54-1540 

33 

33 

Coast Range; Sacramento Valley 

1 , 

3 

CB 

0 

48-1384 

26 

26 

Sacramento Valley; Sierra Nevada foothills 

1, 

3 

CC 

0 

200-3540 

29 

29 

Sierra Nevada foothills and mountains 

2, 

3 




238 

238 





*LAND USE CODE (from: Durrenberger and Johnson. 1976) 

1 Cropland and pasture - Irrigated or uni rri gated 

2 Forest and woodland - grazed or ungrazed 

3 Open shrub woodlano - grazed 

4 Urban 


tested and edited for digitizing and scaling errors. This process 
and the final evaluation will be completed during the next re- 
porting period (June - March 1984). 

1.2.4 Literature Cited 

Daniel. Cuthbert and Fred S. Wood. 1971. Fitting equations to 
data. Wiley- Interscience, New York. 

Durrenberger, Robert W. and Robert B. Johnson. 1976. California - 
patterns on the land. Fifth edition. Mayfield Publishing 
Company. 

1.3 Early Evaluation of TM and MSS Data for Selected Forest Cover Types 
1.3.1 Background 

The objectives for this phase of our research are to: (1) 

develop a basic understanding of the TM data of forested environ- 
ments with regard to spectral characteristics and variability, 
spatial resolution, and radiometric sensitivity, and (2) determine 
the extent to which major forest cover types can be detected and 
Identified on TM digital and Image products. 

Our forestry study site is located In Plumas County, approxi- 
mately 265km northeast of San Francisco. The area contains a 
diversity of forest cover types ranging from stands of red and 
white fir [Abies magnified and A. oonoolor» respectively) to stands 
of mixed conifer dominated by Ponderosa Pine ( Pinue ponderoea ), 
Douglas-flr ( Peeudoteuga memieeii) t and/or Sugar Pine ( Pinue lam- 
bertiana). Several other cover types are also prevalent and Include 
low density Jeffrey Pine (p. jeffreyi) stands on soils derived from 
ultramaflc parent material; hardwood stands; dense shrub fields; wet 
and dry meadows; bare soil; granitic rock outcrops; and large water 
bodies. 


1.3.2 Approach 

The first Landsat-4 scene for this site was acquired on 12 Auoust 
1983 (#84039218143, WRS Path 44, Row 32). The TM data were trans- 
mitted to GSFC via the Tracking and Data Relay Satellite System (TDRSS) 
and the TDRSS receiving station at White Sands, New Mexico. The TM 
data were processed by the Thematic Mapper Image Processing System 
(TIPS) at GSFC as a "P" tape. The simultaneously acquired MSS data 
were purchased from the EDC as a "P" tape (CCT-PM). These data were 
received at U.C. Berkeley on October 14, 1983, on 6250 bpl density 
tapes. 

Large-scale (1:4,300-1:13,000) color and color Infrared obllqjr 
aerial photography was acquired over the study site using a dual 35m>i 
camera system operated from a light aircraft. The purpose of acquiring 
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this photography was to document the environmental conditions 
prevalent at the time of the overpass. These photos are used In 
conjunction with available ground data to document forest canopy, 
non-forest cover, and understory conditions. The aerial photo- 
graphy was acquired within one hour of the Landsat-4 overpass on 
August 12, 1983 from 1200-1300 P.D.T. Ground cover conditions 
were documented based on Interpretation of these photographs, the 
field data sheets from an Intensive field study conducted In the 
area during August 1982, and this Investigator's knowledge of 
forest conditions and management practices In the area. 

Based on the available ground data and aerial photography for 
this site, the TM and corresponding MSS digital data were extrac- 
ted from the full scene CCT's. A 1 200.x 1200 pixel block area 
(1,440,000 pixels) was extracted from the CCT-PT using the IBM 
Image processing system at the Pale Alto Scientific Center with 
assistance of Mr. Ralph Bernstein and Mr. Jeff Lotsplech. The 
corresponding area was extracted from the CCT-PM; and both data 
sets were written to separate disk files for analysis. 

Forest stands of known cover type and condition were located 
In the digital data set, annotated, and statistics calculated for 
each stand. Both radiance numbers (RN), or counts, extracted from 
these sample fields and spectral radiance values were used for sta 
tlstlcal analysis, plots, and Image generation. Spectral radiance 
for each spectral band was calculated based on the following re- 
lationship [1]: 

RN, - 6, 

L, • —4 1 


where, 

•2 -1 -1 

» Spectral radiance, mWcm -sr -um , for TM Band 1 
RN, ■ Radiance number. In counts, for TM Band 1 

A ■ 

0, ■ Average Offset, In counts, for TM Band 1 

* 1 -2 -1 -1 

G, ■ Average Gain, In counts/mWcm -sr xim , for 

1 TM Band 1 
1 ■ Spectral Band \’ndex 

The pre-launch band offset and band gain values were used In the 
equation [2]. The summary statistics for the forestry study site 
are presented In radiance numbers and spectral radiance; forest 
stand statistics and plots are In radiance numbers only. 

The radiance numbers for the entire 1200 x 1200 pixel block area 
were Input to a standard statistical program to compute means, vari- 
ances, coefficients of variation (CV), minimum and maximum values, 
skewness, kurtosls, and covariance and correlation matrices for the 
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YM and MSS data. These summary statistics provide an overview of 
the frequency distributions by band of the TM spectral data, and 
will be used at a later time to optimize color mapping tables to 
produce various Image products for Interpretation. The coefficient 
of variation (CV), expressed as a percent, H used to represent che 
relative variation about the mean, and Is calculated using the mean 
and variance ([std. dev. /mean j x 100). This relative measure allows 
for an easy comparison of spectral band means and their associated 
variability which may exhibit a wide range of values. Tests for 
normality which represent the degree of symmetry and peakedness of 
frequency distributions are skewness and kurtosls, respectively. 

The skewness Is the third product moment about the mean; symmetrical, 
normal frequency distributions have skewness of zero. Large positive 
values of skewness Indicate data non-normality In which an Increasing 
magnitude of the mean Is accompanied by a corresponding Increase In 
variance. Kurtosls Is the fourth product moment about the mean; all 
data distributions have values of kurtosls greater than -2 and normal 
distributions have a kurtosls value of zero. Though both covariance 
and correlation values were calculated, only the correlation matrix 
Is presented. The values of the covariance matrix are difficult to 
Interpret for evaluating the relationship between spectral bands. 

The product moment correlation coefficient (r) Is calculated from 
the band variances, and provides a standard form for evaluating these 
Interband relationships. For the radiance numbers extracted from in- 
dividual forest stands, only means, variances, and CV's were calcu- 
lated. For each TM and MSS spectral band of this site, both black- 
and white and color composite Images were produced to visually repre- 
sent the digital data and to Illustrate the spectral variability of 
selected forest cover types. 

1.3.3 Results and Discussion 

The summary statistics for the forestry study site are shown In 
Tables 2 and 3. The statistics In Table 2 Include the band means, 
variances, CV's minimum and maximum values, and skewness and kurtosls 
for the frequency distribution of the 1,440,000 pixels per band. TM 
Bands 5 and 7 are the most variable bands and can be considered nor- 
mally distributed as Is the least variable band, TM6. TM Bands 1 and 
2 exhibit lower variability but cannot be considered normally dis- 
tributed due to the high values of skewness and kurtosls. No trans- 
formations were performed to normalize the data or stabilize the vari- 
ance. Bands 1, 5, and 6 contain saturated pixels due to (1) the high 
albedo of granitic rock outcrops (Band 1); (2) high albedo of dry, 
bare soil on both level terrain and south facing slopes (Band 5); and 
(3) the high radiant temperature of the granitic rocks and bare soil 
(Band 6). The Interband correlations In Table 3 show the familiar 
high positive correlation between TM Bands 1, 2 and 3; and low corre- 
lation between Band 4 and Bands 1, 2, 3, 5 and 7. Both the short- 
wave bands In the summer forested scene show high positive corre- 
lations with the thermal band. This Is directly related to the close 
association between the radiant temperature response of forest cover 
types, and their relative moisture status. Generally speaking, the 



- 10 - 


ORIGINAL PAGE 13 
OF POOR QUALITY 

TABLE 2. Sunmary statistics for the forestry study site. 


TM 

Band 


_x 

_s 

CV 

Min 

Max 

Skewness 

Kurtosis 

( yJ 

i 

RN: 

66.1 

8.1 

12.3 

41 

255 

3.1 

25.0 


L: 

4.1 

0.5 


2.5 

16.1 



2 

RN: 

23.7 

5.2 

21.9 

11 

166 

2.5 

16.8 


L: 

2.8 

0.5 


1.3 

20.6 



3 

RN: 

23.2 

8.2 

35.3 

7 

214 

2.4 

10.5 


L: 

2.2 

0.8 


0.7 

20.2 



4 

RN: 

64.1 

16.1 

25.1 

0 

172 

0.2 

1.3 


L: 

5.3 

1.4 


-0.06* 

15.8 



5 

RN: 

46.7 

22.5 

48.2 

0 

255 

1.1 

1.5 


L: 

0.6 

0.2 


-0.03* 

3.2 



7 

RN: 

18.0 

11.7 

65.0 

0 

164 

1.5 

3.1 


L: 

0.1 

0.06 


-0.02* 

1.1 



6 

RN: 

190.1 

15.1 

7.9 

139 

255 

1.1 

1.5 


♦noise RN = Radiance number, counts . 

L = Spectral radiance, mWcm~2-sr -um 
n = 1,440,000 pixels 


TABLE 3. Correlation matrix for the seven TW bands for the forestry study 
site. 


BAND 

1 

2 3 

4 

5 

7 

6 

1 

1.00 

0.94- 0.96 

0.15 

0.74 

0.84 

0.30 

2 


1.00 0.96 

0.24 

0.80 

0.87 

0.77 

3 


1.00 

-.18 

0.81 

0.90 

0.81 

4 



,.00 

0.43 

0.18 

0.06 

5 




1.00 

0.93 

0.85 

7 





1.00 

0.97 

6 

n - 

rj,000 pixels 




1.00 
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low reflective properties In TM5 and 7 of vegetative canopies are 
directly related to actively growing vegetation under relatively 
high moisture conditions. As canopy moisture content Increases, 
there Is more absorption of the short-wave radiation by the canopy 
and a correspondingly lower radiant temperature of the canopy due 
to dissipation of heat energy through evapotransplration. Con- 
versely, conditions In a forested environment which have highly 
reflective properties in TM5 and 7 (granitic rock outcrops; dry, 
bare, coarse- textured soil) also have high radiant temperatures 
due to low moisture conditions. 

Individual forest cover type statistics are shown In Table 4 
and plotted In Figure 3. Table 4 contains the band means and 
coefficients of variation (standard deviation/mean x 100) for In- 
dividual vegetated and non-vegetated cover types for the seven TM 
bands. As both sensors were operating simultaneously during this 
overpass, band means and CV's of three of the four MSS bands which 
most closely correspond to the relevant TM bands are also Included 
in Table 4. Though portions of both MSS 3 (0.70-0.81 urn) and MSS 4 
(0.81-1.02 urn) fall within portions of TM 4 (0.78-0.90 urn). MSS 4 
was selected so as to avoid that highly variable region of vege- 
tation response curve between 0.70-0.80 urn where the spectral re- 
flectance of vegetation goes from a minimum at 0.645 un (peak chlor- 
ophyll absorption) to a maximum at 0.880 urn (peak infrared reflec- 
tance). Based on field plot data and synoptic surveys conducted in 
this area for several years, areas of known cover type conditions 
were located in the spectral data and extracted from the disk files 
for statistical analysis. For this phasa of our Investigat’ >n, only 
homogeneous stands of vegetation and other cover types, all of vari- 
able areal extent, were selected for statistical analysis. Future 
analysis will include field data in which specific forejt stand 
parameters (basal area, average stand height, age, diameter breast 
height, timber site) have been measured for 130 plots of variable 
size within the Meadow Valley 7V topographic quadrangle. 

These data from the six reflective bands are graphically displayed 
in Figure 3. For each band, the mean radiance number (in counts), 

+1 standard deviation about the mean, is plotted for each cover type 
sampled. A major cover type discrimination one would like to make 
when stratifying forest cover types is between true fir stands of 
lower economic value and mixed conifer stands of higher economic 
value. By reviewing the spectral statistics, there is no apparent 
spectral difference between these two major cover types. The only 
band in which we see some difference is in the thermal band (TM6). 

This results from the higher topographic position that true fir occu- 
pies. The resulting lower radiant temperature of the true fir canopy 
allows the spectral differentiation between the two types using TM6. 
Another discrimination of importance is between mixed conifer stands 
and stands dominated by hardwood species, namely Black Oak (QueAcua 
ke£Logg<U), which is an indicator of high timber site (highly pro- 
ductive). Statistics from a mixed stand dominated by Black Oak is 
plotted under "Hardwood". From the statistical data, a major dif- 
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FABLE 4. TM and MSS spectral statistics fnr the forestry study site. Values tabulated are average radiance 
numbers (x) and coefficients of variation (CV) calculated from the digital data. 


SPECTRAl 

BAND 

TRUE 

FIR 

x CV 

MIXED 

CONIFER 

x CV 

HARDWOOD 
x CV 

SHRUB 
x CV 

MEADOW 
x CV 

BARE 

SOIL 

X 

CV 

ROCK 
x CV 

WATER 
x CV 

TM1 

63.5 

9.3 

61.6 

3.1 

60.0 

1.7 

65.8 

5.1 

75.1 

10.5 

87.6 

3.1 

99.3 

11.4 

58.3 

3.9 

TM2 

22.2 

14.4 

20.0 

4.8 

; 20.3 

2.8 

25.7 

8.2 

30.6 

14.6 

38.5 


41.8 

14.3 

16.4 

8.4 

MSS1 

15.0 

13.3 

13.4 

9.6 

14.3 

4.0 

16.8 

10.0 

20.7 

13.5 

25.3 


29.0 

10.3 
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Figure 3. Plot of RN values (+ 1 s.d.) for each of the reflective TM bands for 
several forest cover types. 
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ference In spectral response Is seen In TM 4, as one would expect. 

This Is also evident when viewing the in Image products where an 
Increasing percentage of hardwood species within a mixed conifer 
stand Is reflected In Increasing light red tones on the image dis- 
tinguishing these cover types from the "pure" conifer stands having 
a dark red tone. Discrimination of trees from shrubs and herba- 
ceous vegetation Is quite obvious by examining both the spectral 
plots and film products. Differentiation between wet and dry meadows 
can be accomplished using both TM 4 and TM 5 due to the differences 
In high photographic Infrared reflectance (TM 4) and low shortwave 
Infrared reflectance (TM 5) resulting from dense herbaceous growth 
under high canopy moisture conditions In the wet meadows. Dry, 
bare soil and rock outcrops exhibit the highest TM 5 and TM radi- 
ance with some areas on south facing slopes causing pixel saturation 
(RN: 255) on both bands. Deep and relatively clear water bodies in 
this area are spectrally behaving as one would expect with Increasing 
absorption of solar radiation as wavelength Increases. 

These spectral relationships between major cover types are illus- 
trated in Figure 4. Four additive color composites of various TM and 
MSS spectral bands provide a good indication of the relative improve- 
ments of TM data over MSS In terms of both spatial and spectral reso- 
lution. A stereo-pair of small format (35nm), low altitude color 
oblique aerial photography taken coincident with the Landsat-4 over- 
pass Is also shown In Figure 4. The composites cover a ground area 
of approximately 5200 ha (52 km 2 ) and are reproduced here at a scale 
of approximately 1:65,000. Each composite was Imaged directly from 
the digital data using a color graphics camera. Three spectral bands 
were Individually displayed on a high resolution black-and-white moni- 
tor and projected through color filters onto a single frame of 35mm 
Ektachrome (ASA 200) color film. Color Internegatives were then 
generated for making multiple copies of the prints. The top left 
Image (4b) Is a composite of MSS bands 1, 2, and 4; the bottom left 
Image (4c) Is TM 2, 3, and 5; and the bottom right Image (4d) Is 
TM 3, 4, and 6; all projected through blue, green, and red filters, 
respectively. The area coverage of the low oblique aerial photo- 
graphy Is annotated on each composite. Five major forest cover type 
conditions are annotated on the overlay for the upper left TM com- 
posite only. These cover type conditions include: (A) bare soil 

which results from clearing forest vegetation for expanding a pon- 
derosa pine plantation; (B) high density mixed conifer (>70% crown 
closure) Including a variable species mix of Ponderosa pine, Sugar 
pine, White fir, and Incense cedar; (C) low density conifer (Jeff- 
rey pine < 20% crown closure) on soil with low fertility. (D) wet 
(high water table) and dry meadows at/or shortly after peak forage 
production, (E) broadleaf riparian forest vegetation (>90% crown 
closure) Including cottonwood, willow, and big leaf maple. 

The Improvement in the spatial resolution of the TM sensor over 
the MSS sensor can be seen by comparing the top two composites In 
Figure 4 (a,b). The same surface area Is represented on both com- 
posites. Significant detail on the TM Image which cannot be con- 
sistently Identified on the MSS Image Include the location of road 
and stream networks, area and shape of small forest clearings, 
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definition of vegetation boundaries, and discrimination of small 
homogeneous vegetation units within a matrix of complex cover 
types. With the clear definition of major unpaved forest roads 
and discrimination of significant topographical features and 
small stream networks on the TM Imagery, these Image products are 
more suited for field navigation than are the MSS film products. 

A composite which Includes TM 5 (4c) would be best for road and 
stream network delineation because the contrast between a dry road 
surface and moist vegetation surface and between a stream coarse 
and adjacent features Is higher on this band than on a composite 
containing TM 4 without TM 5 (4a). 

The Improvement in spectral resolution of the TM sensor over 
the MSS sensor and the improvements which can be expected by using 
various TM spectral band combinations for stratifying major forest 
conditions can be seen in Figure 4. By comparing the forest cover 
type conditions annotated on the overlay with their appearance on 
the other image types, several significant results are evident. 

Image composites which contain TM Band 5 In combination with two 
other visible bands (i.e. TM 2, 3) are inferior than those which 
contain TM 4 for the discrimination of important highly infrared 
reflective (TM 4) cover types (shrubs, hardwoods, meadows). This 
is evident by examining the detail in Figure 4a (TM 2, 3, 4) in 
which a variety of these cover types of varying areal extent are 
easily dlscriminatea due to the bright red tones. By examining 
the same detail In Figure 4c (TM 2, 3, 5) these cover types are 
not as evident. This includes the riparian for*»st area at "E", a 
network of small wet meadows near the upper left center portion of 
the image, and the dense shrub field in the upper rig.tt center of 
the image. This results from the high level of absorption of the 
solar radiation in this band (1.57-1.78 urn) by vegetation having 
a high leaf water content. The optimum combination of reflective 
TM bands for discriminating most forest cover types would be either 
a TM 2, 4, 5 or TM 3, 4, 5 where both the reflective and absorptive 
properties of diverse cover types could be exploited. 

The discrimination of thermal differences of forested cover 
types is important for assessing regional soil and plant moisture 
conditions and for mapping those sites suitable for reforestation 
programs. The data being acquired by the broad thermal infrared 
band, TM 6 (10.42-11.66 urn), is the highest quality data both 
spatially and spectrally, ever acquired by an earth-orbiting land 
remote sensing satellite. The TM 6 data have added a new dimension 
to the characterization of forest resources on image products. This 
is evident on the composite shown in Figure 4d. This composite in- 
cludes two reflective bands (TM 3, 4) and the thermal band (TM 6). 
Though the coarser ground resolution (120m vs. 30m) prohibits the 
discrimination of small features of interest (small forest clearings, 
water bodies), the large scale radiant temperature differences between 
major forest cover types is obvious. In this composite, those cover 
types which are highly infrared reflective in TM 4 and also have a 
high moisture status (low radiant temperature) are shown in tones of 
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light green (wet meadows, riparian forests, dense shrub fields). 

At the other extreme, those cover types with relatively low In- 
frared reflectance In TM 4 and also have a low moisture status 
(high radiant temperature) are shown In tones of light red, or 
pink (bare soil, dry meadows, low density coniferous forest). In- 
termediated cover types, such as high density conifer forest on 
northerly aspects, are shown In darker green tones. By combining 
the thermal band of low spatial resolution with the reflective 
bands of high spatial resolution, a unique and useful Image product 
Is produced. This "thermal" composite should Include TM 4 as one 
of the reflective bands and not TM 5. Due to the high positive 
correlation between TM 5 and^TR 6 (Table 3.) resulting from the 
radiant temperature response and the moisture relationships of 
diverse cover types, the composite containing both TM 5 and TM 6 
would be redundant and provide less Information about forest cover 
conditions than would a composite containing the TM 4 band. Again, 
the most useful composite would be one containing both TM 4 and TM 5, 
or both TM 4 and TM 6 with the addition of a visible band (TM 1, 2, 3). 

1.3.4 

Spectral data were simultaneously acquired by the Thematic Mapper 
(TM) and Multispectral Scanner (MSS) aboard the Landsat-4 spacecraft 
on August 12, 1983 (Path 44) during an operational TDRSS test. The 
data received from TDRSS were processed through the Thematic Mapper 
Image Processing System (TIPS) at NASA-Goddard Space Flight Center. 

A systematic analysis of both image and numeric data for a forestry 
study site In the northern Sierra Nevada, California, has yielded 
the following significant results: 

(1) TM Bands 1, 5, and 6 contain saturated pixels (RN = 255) 
for areas of high albedo and high radiant temperature (bare 
soil, rock outcrops). 

(2) TM Bands 5, 7, and 6 show high positive Interband correla- 
tion due to the close association between the absorption of solar 
radiation in the 1.57-1.78 urn and 2.10-2.35 urn wavelength bands 
by vegetation and the low radiant temperature of the same vege- 
tation due to the cooling effects of evapotransplration. 

(3) The TM reflective bands (TM 1-5, 7) do not permit a clear 
discrimination between true fir and mixed conifer stands on 
spectral reflectance characteristics alone. Use of both the 
thermal band (TM 6) and elevatlonal gradients will improve this 
Important discrimination. 

(4) The high spectral, spatial, and radiometric resolution of 
the reflective TM bands permit the discrimination of forest 
stands with Increasing densities of Black Oak, an Indicator 

of highly productive forest sites. 
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(5) The Improved spatial resolution of the TM sensor over 
that of the MSS allows for Improved discrimination of road 
and stream networks, area and shape of small forest clearings, 
significant vegetation boundaries, and small homogeneous vege- 
tation units within a matrix of complex cover types. 

(6) Color Image composites which Include TM 5 are best suited 
for discriminating these small forest features essential for 
both navigation and mapping purposes. 

(7) Color image composites which Include TM 5 with other 
visible bands (l.e. TM 2, 3) are Inferior for discriminating 
Important vegetation types which are highly reflective In the 
photographic infrared portion of the spectrum (TM 4, 0.78-0.90 urn). 

(8) The optimum combination of reflective TM bands for discrim- 
inating most forest cover types would be either TM 2, 4, 5 or 

TM 3, 4, 5 where both the reflective and absorptive properties 
of diverse cover types can be exploited. 

(9) The thermal Infrared band (TM 6) contains valuable Infor- 
mation on the relative radiant temperature differences between 
major cover types. Image products combining the spectral 
qualities of both the thermal and reflective bands and the 
spatial qualities of the reflective bands provide an optimum 
multispectral approach for the stratification of major forest 
cover type conditions. 

(10) Color image composites should not contain both TM 5 and 
TM 6 due to the high positive correlation between these bands 
in a forested environment. If spatial resolution must be maxi- 
mized, a composite containing TM 5 Instead of TM 6 Is recom- 
mended. 


In summary, the spectral, spatial, and radiometric quality of 
the TM data are excellent. Based on this early stage of research, 
we conclude that this high quality data are more than sufficient for 
meeting most of the Inventory objectives of the renewable resource 
specialist. The TM data should prove extremely valuable for (1) estl 
mating forest cover types and area proportions required for most In- 
ventory and sampling purposes; (2) updating land use survey maps, and 
(3) determining the size, shape, and location of Individual forest 
clearings and water sources. 
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Figure 4. Landsat-4 TM and MSS color (additive) composites and color aerial oblique photography 
illustrating major forest cover type conditions. 4a: TM Bands 2,3,4; 4b: MSS Bands 
1,2,4; 4c: TM 2,3,5; 4d: TM 3,4,6. The bottom image (4e) is a stereogram of a portion 
of area covered by the composites using the low altitude (600m above the terrain) 35mm 
color aerial oblique photography acquired coincident with the Landsat overpass on 
August 12, 1983. The polygon on the composite overlay depicts the area covered by the 
two oblique photographs. See the text for the cover type conditions annotated on the 
overlay. 
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Brochure. 

2.2 Presentations 

PECORA VIII Memorial Symposium, Sioux Falls, South Dakota, 

October 4-7 (S. DeGloria) 

California Map Society Meeting, November, San Francisco 
(A. Benson) 

Remote Sensing Training Course, Mendocino National Forest, 
Willows, California, November 8-9 (A. Benson, K. Durnner, 

S. DeGloria) 

Third Landsat-4 Investigator's Workshop, NASA-Goddard Space 

Flight Center, Greenbelt, Maryland, December 6-7 (S. DeGloria) 

Comnunity College Instructor's Workshop, Modesto, California, 
December 10 (B. Wood) 

Distinguished Visitor's Seminars (S. DeGloria): 

Mr. Camino, Director General, Promotion of Agriculture 
for Spain 

Mr. Nedal, Director, Spanish Agricultural Research Service 

Mr. Tasias, Directjr, Spanish Agricultural Experimental 
Station at Tarragona 

Mr. Mathura D. Rajbhandari, Dean, Institute of Forestry, 
Hetauda, and Tribhuwan University, Tritureswor, Katmandu, 
Nepal 


3.0 FUNDS EXPENDED TO DATE 


The funds expended to December 31, 1983 under this contract are summarized 
in NASA Form 533M, "Monthly Contractor Financial Management Report", dated 
January 31 , 1984. 


4.0 PROBLEMS ENCOUNTERED TO DATE 

Specific problem areas for this reporting period are discussed in the 
October, November and December monthly reports submitted to Mr. Darrel Williams, 
Code 923, NASA-Goddard Space Flight Center. 
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